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Abstract 

Oscillations of magnetic structures in the solar corona have often been interpreted 
in terms of magnetohydrodynamic waves. We study the adiabatic magnetoacoustic 
modes of a prominence plasma slab with a uniform longitudinal magnetic field, sur- 
rounded by a prominence-corona transition region (PCTR) and a coronal medium. 
Considering linear small-amplitude oscillations, the dispersion relation for the mag- 
netoacoustic slow and fast modes is deduced assuming evanescent-like perturbations 
in the coronal medium. In the system without PCTR, a classification of the oscilla- 
tory modes according to the polarisation of their eigenfunctions is made in order to 
distinguish modes with fast-like or slow-like properties. Internal and external slow 
modes are governed by the prominence and coronal properties respectively, and fast 
modes are mostly dominated by prominence conditions for the observed wavelengths. 
In addition, the inclusion of an isothermal PCTR does not substantially infiuence 
the mode frequencies, but new solutions (PCTR slow modes) are present. 
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1. Introduction 

Despite recent observational and theoretical advances (see the reviews by Martin, 
1998; Patsourakos and Vial, 2002), the intimate physical nature of solar prominences 
is still enigmatic. It is clear that prominences can only exist at locations where the 
coronal magnetic field is "peculiar" enough to allow for the equilibrium of the dense 
prominence plasma against the force of gravity. Nevertheless, there remain many 
unknowns in the birth, evolution, and death of these objects. 

Prominence seismology attempts to contribute to the research of prominence 
properties. Observations of waves and oscillations in prominences have become in- 
creasingly complex in the last decade and a rich field has emerged. Oscillatory 
periods ranging from less than a minute (Balthasar et al, 1993) to 12 hours (FouUon, 
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Verwichte, and Nakariakov, 2004) have been detected. It has also been found that 
perturbations are typicaUy attenuated with a damping time of a few periods (Ter- 
radas et al, 2002) and in some cases the phase speed and wavelength of propagating 
perturbations have also been characterised {e.g. Terradas et al, 2002). Finally, al- 
though there is a considerable gap between observations and theory, some attempts 
of doing prominence seismology have been undertaken and reasonable values of 
the prominence physical parameters have been obtained (Pouget, Bocchialini, and 
Solomon, 2006). The observational background has been reviewed at length in various 
publications (Oliver and Ballester, 2002; Wiehr, 2004; Engvold, 2004; Banerjee et al., 
2007) to which the reader is referred for more information. 

The theoretical magnetohydrodynamic models developed so far start from very 
crude representations of the prominence and usually focus on a few "ingredients" 
whose relevance wants to be tested. For example, Joarder and Roberts (1993) studied 
the adiabatic modes of oscillation of a plasma slab threaded by a skewed magnetic 
field. In Oliver and Ballester (1996) the magnetic geometry was simplified by taking 
a transverse magnetic field to the prominence slab, but a prominence-corona transi- 
tion region (PCTR) was included; this is the only theoretical study of prominence 
oscillations in which the infiuence of the PCTR has been taken into account. The 
temporal and spatial damping of perturbations has been investigated in some other 
works [e.g. Terradas et al, 2005). Oliver and Ballester (2002) and Banerjee et al. 
(2007) kave reviewed this subject. 

The purpose of this work is to address the importance of the PCTR on the 
oscillations of a prominence slab threaded by a straight, uniform magnetic field 
parallel to the slab axis (and not perpendicular to it as in Oliver and Ballester, 1996). 
The PCTR is a narrow layer through which the physical parameters (temperature, 
density, ionisation degree, etc.) change abruptly from prominence to coronal values 
(Cirigliano, Vial, and Rovira, 2004). The availability of telescopes for the observation 
of spectral lines sensitive to PCTR temperatures has allowed the detection of waves 
in this region (Blanco et al, 1999; Bocchialini et al, 2001; Pouget, Bocchialini, and 
Solomon, 2006). The first of these works reports the propagation of a disturbance 
travelling at w 170 km s~^ and with a period in the band 3.2-6.4 minutes. Pouget, 
Bocchialini, and Solomon (2006) reported periodicities detected in three solar fila- 
ments in the He I line at 5843.3 A (corresponding to 20000 K) and identified some 
of the periods with the six fundamental oscillatory modes according to the model 
proposed by Joarder and Roberts (1993). Incorporating the PCTR in our model will 
allow us to test whether it has an observable effect on the oscillations of a prominence. 

When the PCTR is removed, our equilibrium configuration is identical to that 
of Edwin and Roberts (1982), hereafter ER82. ER82 studied a wide range of values 
for the physical conditions of the slab and the outer medium, plotting dispersion 
diagrams for each situation, and also discussed the dispersion relation in the long- 
and short-wavelength limits. Joarder and Roberts (1992), hereafter JR92, extended 
the work of ER82 by considering the possibility of motions and propagation in the 
2/-direction and made a more refined classification of the oscillation modes according 
to their phase speed. Our purpose is to revise the work done by JR92 in the case of 
no motions and no propagation in the y-direction (i.e. the case analysed by ER82), 
taking into account realistic physical conditions similar to those in a prominence and 
its surrounding PCTR and corona, and to study the polarisation of motions and the 
other perturbations. Additionally, we explore the parameter space in an attempt to 
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Figure 1. Sketch of the equilibrium. 

understand how the variation of tlie temperature, density and magnetic field afi^ects 
the oscillatory period. 

This paper is organised as follows: Section 2 contains a description of the equi- 
librium model and the linear, adiabatic wave equations, besides the derivation of 
the dispersion relation of magnetoacoustic modes. The oscillatory modes are studied 
first for a system without PCTR (as a revision of ER82) in Section 3. Next, an 
isothermal PCTR is included in the equilibrium and its effects are investigated in 
Section 4. Finally, our conclusions are given in Section 5. 



2. Equilibrium, Basic Equations, and Dispersion Relation 

Our equilibrium configuration (Figure 1) is made of an homogeneous plasma layer 
with prominence conditions embedded in an unbounded corona with a prominence- 
corona transition region between them. The PCTR is represented as an isothermal 
and homogeneous layer. The physical conditions (two different temperatures for 
the PCTR have been considered) together with the characteristic velocities of each 
medium are displayed in Table 1. PCTR and coronal densities are computed by fixing 
their respective temperatures and imposing pressure continuity across the interfaces. 
The magnetic field is orientated along the a;-direction, having a value of 5 G in all 
layers. The three media are unlimited in the x- and y-directions. The half-width of 
the prominence slab is z-p = 3000 km and the PCTR width is Atr = 900 km. 

In our study of linear, adiabatic, magnetohydrodynamic waves we follow the 
derivation of ER82 (see also Roberts, 1981) and use their Equation (5), which governs 
the propagation of magnetoacoustic fast and slow modes. The Cartesian axes have 
a different orientation in ER82 and here, so ER82's Equation (5) is now written 



adiabatic.tex; 2/02/2008; 4:19; p. 3 



R. Solor, R. Oliver, and J.L. Ballcster 



with 



-——+kzVz=0, (1 



where Vz is the normal velocity component, uj is the frequency, kx and kz are the 
components of the wavemimber, Cg = is the adiabatic sound speed squared, 

the Alfven speed squared and c-p is the cusp (or tube) speed squared, 

= (3) 

The parameters Bq, pq, and po represent the equilibrium magnetic field strength, 
pressure, and density of any region. In all the following formulae in this work, sub- 
scripts "p", "tr", or "c" denote quantities computed using prominence, PCTR, or 
coronal values, respectively. 

The expressions for the perturbed quantities as a function of Vz and its derivative 

are 

(4) 



ikxCs 


dvz 




dz ' 


iujpo 


dvz 


uJ^ — kxci 


dz ' 


iupocs 


dvz 



~ c^-kld^' 

^ it^PQCs dvz 

u'^ — kxci dz 

_ iBodvz , 

Bix -7—, [") 

ui dz 

Biz = Vz, (8) 

with Vx the velocity perturbation along the magnetic field and pi, pi, Bix, and 
B\z the perturbations to the density, pressure and the x- and z-components of the 
magnetic field. Next, we write expressions for the perturbations to the magnetic 
pressure (pim) and the total pressure (pit), 

_ Bo _ ipQv'i dvz I s 

Plm — Bix — (9) 

p, UJ dz 

ipo {u^ - klvl) dVz 

Pit = Pi+Pim = ^ 72 --7-- (10) 

ujki dz 

In order to obtain an analytical dispersion relation for the linear adiabatic magne- 
toacoustic waves, we must impose a solution that satisfies Equations (1) and (2) for 
our equilibrium. Since we are interested in wave modes which are evanescent in the 
corona we impose the restriction k^^ < 0. Regarding the prominence, solutions with 
kzp > correspond to body waves, whereas solutions with fe^p < are surface-like 
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Table 1. Temperature, density, and mean atomic weight together with the 
cusp, sound, and Alfven speeds for all the considered media. Velocities are 
expressed in kms~^. 
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pa (kg m 
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CT 


Cs 




Prom. 


8 000 


5.00 X 10" 


11 


0.8 


11.56 


11.76 


63.08 


PCTR (1) 


50 000 


5.00 X 10- 


12 


0.5 


36.56 


37.19 


199.47 


PCTR (2) 


500 000 


5.00 X 10" 


13 


0.5 


115.62 


117.62 


630.78 


Corona 


1 000 000 


2.50 X 10" 


13 


0.5 


163.51 


166.33 


892.06 



modes. On the other hand, in the transition region the case fc^tr > corresponds 
to a propagating wave, whereas the case fc^tr < corresponds to an evanescent, 
tunneUing wave in tlie PCTR. All of these cases are allowed and studied in this 
investigation. The analytical function proposed for Vz as solution of Equation (1) is 
then a piece-wise function, 



f —2c < z < — ^p, 
f -Zp < 2 < 2p, 
f 2p < 2 < 2c, 



Ai exp [fc^c (2 + 2c)] , if 2 < -2c, 

A2 cos [fcztr (2 + 2p)] + A-i sin [fcztr (2 + 2p)] , 

Vz{z) = ^ Ai cos {kzp 2) + Ar^ sin (fc^p 2) , 

cos [fcztr {z - 2p)] + A7 sin [fcztr {z - 2p)] , 

An exp [-kzc (2 - Zc)] , if 2c < 2. 

(11) 

This solution satisfies the evanescent assumption, — > when 2 — > ±cxd. Imposing 
continuity of Vz and the total pressure perturbation across the interfaces, we obtain 
eight homogeneous linear algebraic equations for the constants Ai — Ag. The non- 
trivial (i.e. non-zero) solution of this system gives us the dispersion relation 

^ (klvj^tr - kzp I I (fczp2p) TZi ± (klv\p - Uj'^^ fcztr 7^2 = 0, (12) 



with 



n 



PckztT (klv\c - u'^^ cos [fcztr [Zc - 2p)] + 
Ptrkzc (klv\t^ - uj'^^ sin [fcztr {Zc - 2p)] , (13) 
1^2 = ptikzc (^klv\ti- - t^^) COS [fcztr (2c - 2p)] - 

Pckzti- (klvAc - sin [fcztr (zc - Zp)] , (14) 

where cot/tan terms and ± signs in Equation (12) are related with the symmetry of 
the perturbations: the cot term and -|- sign correspond to kink modes {A5 = 0), while 
the tan term and — sign correspond to sausage modes [A^ = 0). When the PCTR 
is removed from the equilibrium (2c = 2p), the dispersion relation (12) reduces to 

y- (fc^fic - kzp I I {kzpZp) ± (^klv\p - uj'^^ kzc = 0. (15) 

This dispersion relation is equivalent to Equation (11) in ER82. 
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3. Results: Equilibrium without PCTR 

First, we consider an equilibrium without PCTR and take ky = 0, which is equivalent 
to that of ER82. The study of this simplified system allows us to understand the 
basic solutions of the dispersion relation and their behaviour, which will be useful 
in order to compare with the results obtained when a PCTR is included in the 
system. Since ER82 and JR92 investigated in detail this equilibrium, we concentrate 
our attention in the study of the oscillatory modes according to the polarisation of 
their eigenfunctions and to their response when the equilibrium physical parameters 
(temperature, density, magnetic field, etc.) are changed. 

3.1. Phase Speed Diagram 




10'" 10'^ 10'' 10'' 10° 10' 10' 10' 



Figure 2. Phase speed (in km s~^) versus the dimensionless longitudinal wavenumbcr of the 
normal modes of a prominence-corona system without transition region. Solid lines denote 
kink modes and dotted lines denote sausage modes. Only the fundamental modes are plotted. 
The shaded zones are forbidden regions where evanescent-like solutions in the corona do not 
exist. Note that the vertical axis is not drawn to scale. 

In Figure 2 the magnetoacoustic wave dispersion diagram is displayed. As noted 
in ER82, the present ordering of Alfven and sound speeds implies that externally 
evanescent solutions only exist in three phase speed ranges (or windows): cxp < 
u)/kx < Cgp, WAp < i^/kx < cxc and Csc < uj/kx < wac- Each window contains 
an infinite number of harmonics, although we restrict our study to the fundamental 
kink and sausage modes and so only two curves are represented in each window 
of Figure 2. In the case ky = 0, JR92 denoted that the modes in the three phase 
speed ranges as slow, fast externally slow and fast externally fast (or fast ES and 
fast EF for short), respectively. It is worth mentioning that for ky = and for 
the present ordering of Alfven and sound speeds this equilibrium does not support 
surface-like modes {kip < 0) since all solutions are of the form of body waves (/c^p ^ 
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0). Notice that both ER82 and JR92 missed the horizontal (i.e. non-dispersive) 
solutions in the top two windows of the phase speed diagram (see their Figure 7 and 
Figure 2(a), respectively). The study of the nature of the solutions carried out in the 
next two subsections suggests that the classification of JR92 for ky = might not 
fuUy reflect the physical nature of the modes. A more representative description of 
the magnetoacoustic modes is as follows: 

• Solutions in the bottom window are internal slow modes, that exist for all values 
of the longitudinal wavenumber (kx) and whose features are mainly determined 
by the prominence physical properties. The phase speed of the internal slow 
modes is enclosed in a very narrow region between cxp and Cgp, so a good 
approximation for the frequency is ~ Cgp kx • 

• The dispersive solutions in the upper two windows are the kink and sausage fast 
modes. The phase speed of fast modes varies between v\p and t;ac i but they are 
separated in two different windows by the "forbidden region" (ctc < uj/kx < 
Csc), which causes the presence of lower and upper cut-off frequencies. In this 
forbidden region, evanescent waves are not possible since solutions become leaky 
{kzc > 0), which means that they freely propagate in the coronal medium and 
carry energy away from the prominence slab. The "forbidden region" causes the 
appearance of an horizontal branch with uj/kx ^ Ctc, which corresponds to a 
solution whose properties are slow-like and fixed by the coronal medium, as we 
explain later in Section 3.3. For this reason, we call this branch of the solution 
the external slow kink mode, although one must bear in mind that this mode is 
not propagating in the corona and the word "external" only points out that its 
behaviour is dominated by the properties of the outer medium. In Figure 3 we 
show an enlargement of the zone in the dispersion diagram close the forbidden 
region, where the fundamental fast kink mode separates into two branches and 
the external slow kink mode arises. 
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In addition, Figures. 2 and 3 show the presence of a cut-ofi frequency in tire upper 
branch of tire external slow kink mode: this solution is only evanescent for kxZp < 0.6 
and becomes leaky (fc^^, > 0) for kxZp > 0.6. Therefore, it is not drawn in Figures. 2 
and 3 for kxZp > 0.6. The cut-off frequency is present in the upper branch of the 
external slow mode and not in the lower branch due to the fact that its phase speed 
tends to cluster towards the external tube speed, ctc- The lower branch approaches 
ctc as an evanescent-like solution since the leaky region is not crossed. On the other 
hand, the upper branch of the external slow solution becomes leaky because the sign 
of fc^c changes from negative to positive when its phase speed reaches the external 
sound speed (see Equation (2)) in its way to the clustering towards ctc- 

3.2. Properties of the Oscillatory Modes 

In order to study in depth the fundamental oscillatory modes and their magne- 
toacoustic properties, we plot the eigenfunctions for some selected values of the 
wavenumber. In particular, we pay attention to the perturbed velocities {vz and 
Vx)j the perturbed density (pi), and the x-component of the perturbed magnetic 
field (Bix) (see Figure 4). The behaviour of non-leaky kink and sausage solutions is 
similar, so we concentrate on the former. 

Beginning with internal slow modes, they do not have cut-off frequencies and so 
they exist for all values of kxZp. The eigenfunctions corresponding to the fundamental 
slow kink mode are plotted only for kxZp = 0.1 in Figure 4 (top row) since they are 
similar for other values of kxZp- We can observe that perturbations are essentially 
confined to the slab, and the dominant component of the perturbed velocity is the 
one parallel to the magnetic field {vx)- On the other hand, pi and Bix are also large 
in the slab, as for slow waves. Since perturbations are much larger in the prominence 
than in the corona, we expect that internal slow modes are dominated by prominence 
conditions whereas the corona has little or no influence. 

Regarding fast modes, they are dispersive and their phase speed varies between 
WAp and v\c- All of these modes, except for the fundamental kink, have an upper cut- 
off frequency (ojcut)- For uj > ci.'cut, modes become leaky. Moreover, additional cut-off 
frequencies occur due to the presence of the forbidden region cxc < c^/kx < Csc- The 
eigenfunctions of the fundamental fast kink mode for kxZp = 0.01 and kxZp = 1 can 
be compared in Figure 4 (second and third rows). The confinement of perturbations 
to the prominence slab is poorer than for the slow modes and decreases with de- 
creasing kxZp. The perturbations for small kxZp have long tails which penetrate far 
into the corona, so we expect that these solutions are highly infiuenced by coronal 
conditions. This is in agreement with Diaz et al. (2001), who studied fast oscillations 
in a Cartesian line-tied slab and found that the amplitude of transverse oscillations 
can still be large in the corona at a distance of many slab widths. Such as corresponds 
to a fast mode, the dominant component of the velocity is Vz and the density and 
longitudinal magnetic field perturbations are much smaller than in the case of slow 
modes. For kxZp = 1, both components of the velocity have a similar amplitude in 
the corona because the fast mode is externally slow. However, Vz is dominant in the 
corona for kxZp = 0.01 because the fast mode is now externally fast. 

Finally, the perturbations for the external slow mode are plotted for kxZp = 0.01 
and kxZp = 0.5. This solution presents a strong, slow-like polarisation outside the 
slab, with the amplitude of Vx much larger than that of Vz, and produces significant 
density and longitudinal magnetic field perturbations compared to those of the fast 
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Figure 4. Perturbations v^,, Vx, pi, and Bix (in arbitrary units) versus the dimensionless 
distance to the prominence axis corresponding to the fundamental kink modes for some selected 
values of kxZp. The shaded region corresponds to the prominence slab. The plots in the second 
row (fast mode for kxZp = 0.01) have a different scale for the horizontal axis to emphasize the 
poor confinement of this solution. 





mode. Nevertheless, the polarisation of oscillations in the prominence slab is fast- 
like, i.e. Vz is larger, but since the amplitude of perturbations is much larger in the 
corona than in the slab, one expects that coronal properties govern the behaviour of 
this solution, its slow-like character being therefore dominant. This is verified in the 
next Section. 
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3.3. Dependence on the Equilibrium Physical Conditions 

Next, we investigate how the phase speed of the fundamental modes is affected by 
changing the equilibrium physical parameters. For simplicity, we limit this analysis 
to kink modes (Figure 5). 
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Figure 5. Phase speed of the fundamental slow kink mode with fc^Zp = 0.1 (solid line in the 
top row) and the fundamental fast and external slow kink modes with kxZp = 0.1 (dotted and 
dashed lines in the bottom row, respectively). From left to right the four columns display the 
variation of the phase speed against the prominence temperature, the prominence density, the 
magnetic field and the coronal temperature. 



In the case of the internal slow mode for kxZ-p = 0.1, the phase speed only depends 
on the prominence physical conditions whereas the value of the coronal temperature 
does not affect it. The dependence of the phase speed on the temperature can be 
understood from the above mentioned approximation for the internal slow mode 

1 /2 

frequency iuj « kxC^p), leading to tiph ~ Csp, so Wph ~ 7p . Finally, we see that 
the phase speed of slow modes decreases for weak field and large density. This fact 
can be understood taking into account that the phase speed diagram presented in 
Figure 2 is modified for weak magnetic field and large density. In this situation, ii^p 
becomes smaller than Csp and then the windows of existence of non-leaky fast and 
internal slow modes overlap, and both modes couple. This causes the behaviour of 
the internal slow mode phase speed for weak field and large density. 

Regarding the fast and the external slow modes, their phase speeds for kxZ-p = 0.1 
are plotted together in Figure 5 (bottom panels). The fast mode is almost indepen- 
dent of prominence and coronal temperatures, and the weak dependence that is 
shown is due to the closeness of the coupling with the external slow mode. The 
values of the prominence density and magnetic field affect the fast mode in an 
important way. The phase speed increases with the magnetic field in a linear way, 
whereas it decreases when large values of the prominence density are considered. Both 
dependencies are consistent with the rough approximation iipjj ^ 'L'Api so Vp\^ ~ Bq 

and Wph ^ Pp 

Finally, the external slow mode only depends on the temperatures, but it is 
unaffected by the value of the prominence density and magnetic field. The phase 
speed slightly decreases when the prominence temperature grows, whereas it grows 
dramatically with the coronal temperature according to the approximate relation 
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iipij ~ Tc , SO one can conclude that the corona is the dominant medium regarding 
tlie properties of this mode, as we suggested in Section 3.2. In fact, the external slow 
mode behaves as the slow mode in an unbounded medium with coronal properties, 
although we must bear in mind that it is not a propagating but an evanescent wave 
in the external medium. In addition, this solution does not exist for small density 
and strong field because the external tube speed, cxc, is then less than ti^p under 
these conditions. 



4. Results: Equilibrium with an Isothermal PCTR 

4.1. Modification of the Phase Speed Diagram 

Now, we include a PCTR in the ecjuilibrium configuration. Our aim is to study 
the effects arising from the presence of the PCTR on the basic oscillatory modes 
of the system without PCTR, which have been investigated and classified in the 
last section. Since we assume an isothermal transition between the prominence and 
coronal media, the PCTR temperature (Ttr) is a relevant parameter. One must bear 
in mind that to consider an isothermal PCTR is, obviously, a rough approximation 
but it allows us to obtain an analytical expression for the dispersion relation. 

When a PCTR is present in the system, solutions are propagating in the PCTR 
{kfti > 0) if their phase speed lies in the region cxtr < i^/kx < Csu-. Outside this 
region, solutions are tunnelling in the PCTR {k'^tr < 0). Depending on the value of 
the PCTR temperature, the character of the solutions can be different, such as is 
summarised next. 

• If Cstr < VAp (Ttr < 144000 K for our equilibrium parameters), a new set 
of solutions of the dispersion relation (Equation (12)) is found in the region 
CTtr < ijj/kx < Cstr, which falls into the zone Csp < uj/kx < VAp, where the 
solutions are surface-like in the prominence (fc^p < 0). This is the only case in 
which surface waves are possible in the prominence for ky = 0. 

• On the contrary, for cxtr > I'Ap (Ttr > 149 000 K) the region cxtr <oj/kx < Cgtr 
is embedded into the middle propagating window in the dispersion diagram, 
and so solutions with cxtr < (^/kx < Cgtr are now body-like (fc^p > 0) in the 
prominence. 

For simplicity, we define T^^ as the PCTR temperature which separates one situation 
from the other. It is worth mentioning that an intermediate situation exists when 
CTtr < VAp < Cstr and although it could be of academic interest, it corresponds to 
a very peculiar case and is not examined here. Since taking a uniform Ttr in the 
equilibrium is an idealisation, both Ttr < Tt* and Ttr > Tj* are equally (un)realistic 
and so both must be considered. Nevertheless, the presence ofaPCTRwithTtr <Tt*r 
does not modify in an important way the phase speed of fast and slow modes obtained 
for the system without transition region (compare Figure 6(a) with Figure 2) since 
only a new isolated window of solutions (cxtr < i^/kx < Cgtr) is added to the phase 
speed diagram. This new window contains now a new kind of wave modes which 
we label PCTR slow modes, since they show a clear slow-like character and are 
dominated by the PCTR physical conditions (see Section 4.2). An infinite mmiber of 
kink and sausage harmonics is enclosed in this thin window. These new solutions are 
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Figure 6. Phase speed (in km s~^) versus the dimensionless longitudinal wavenumbcr of the 
normal modes of a prominence-corona system with transition region. Solid lines denote kink 
modes and dotted lines denote sausage modes. Only the fundamental modes are plotted. The 
shaded zones are forbidden regions where evanescent-like solutions of the dispersion relation do 
not exist. Computations performed for a PCTR temperature of a) 50 000 K and b) 500 000 K. 
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non-dispersive propagating waves in the transition region {k^tr > 0) but surface-like 
in the prominence (fc^p < 0) 

On the other hand, the situation Ttr > Tj* has more physical consequences and, 
therefore, is investigated in more detail next (see Figure 6(b)). For Ttr > Tj* the 
PCTR slow modes contained in the region cxtr < L^/kx < Cstr interact with the 
fast modes by means of couplings. In this situation, the PCTR slow fundamental 
kink mode has an special relevance because its behaviour for small kxZp is different 
from that of the other PCTR slow modes. To illustrate this especial fact, a zoom 
to the dispersion diagram is performed in Figure 7, displaying only the fundamental 
kink modes. There, we see that the fast mode and the PCTR slow mode show 
a coupling for kxZp w 0.4. At the coupling, the curves show an avoided crossing 
and the oscillatory modes swap their magnetoacoustic properties. On the left-hand 
side of the coupling, the PCTR slow mode is a tunnelling solution in the PCTR 
{kttj: < 0)i whereas on the right-hand side this solution is body-like {k'^tr > 0)- This 
peculiar behaviour of the PCTR slow fundamental kink mode is not observed for the 
other PCTR slow modes, which are always propagating waves enclosed in the region 
CTtr < ijj/kx < Cstr. The behaviour of the PCTR slow harmonics can be seen in more 
detail in Figure 8. The phase speed of these PCTR slow harmonics clusters towards 
CTtr 5 while the fast modes are dispersive and couple with the PCTR slow harmonics 
when their frequencies coincide. With respect to the prominence, all these solutions 
are body-hke (fc^p > 0). 




Figure 7. Zoom of Figure 6(b) to show the couplings of the fundamental fast kink mode (solid 
line) with the external slow kink mode (solid line and O) and the PCTR slow kink mode (solid 
line and A). Note that the PCTR slow kink mode has a propagating branch (fc^jj. > 0) inside 
the region cxtr < '^/kx < Cgtr (light-shaded region), and a tunnelling or evanescent branch 
i^ltr < 0)- The dark-shaded region (cxc < (^/kx < Cbc) corresponds to the forbidden, leaky 
band. The phase speed is given in km s^^. 

In this way, we notice some parallelism in the behaviour of the PCTR slow fun- 
damental mode and the external slow mode, since their phase speeds tend to cluster 
towards the tube speed of their respective dominant medium (cxc for the external 
mode and cxtr for the PCTR mode). Nevertheless, there is an important difference 
between their behaviours because the external slow mode cannot penetrate as a non- 
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leaky solution in the region cxc < i-o/kx < Csc, and so its existence is restricted to 
small values of kxZp. 




Figure 8. Behaviour of the kink-hke solutions inside the narrow region cxtr < t^/kx < Cetr 
(shaded zone) for a PCTR temperature of 500 000 K. The quasi-horizontal lines correspond to 
PCTR slow modes whereas the dispersive lines are fast modes. Only several low harmonics are 
plotted, since the shaded region contains an infinite number of PCTR slow harmonics whose 
phase speeds cluster towards cxtr- The phase speed is given in km s~^. 



4.2. Magnetoacoustic Properties of the New Oscillatory Modes 

As we have just described, the presence of an isothermal PCTR produces a mod- 
ification of the phase speed diagram and new solutions, the PCTR slow modes, 
are present. To check their magnetoacoustic character, we fix the PCTR tempera- 
ture to 500 000 K (case Ttr > T^*) and plot the eigenfunctions of the PCTR slow 
fundamental kink mode in Figure 9 for kxZp = 0.25 (tunnelling branch) and for 
kxZp = 1 (propagating branch). Perturbations are essentially confined within the 
PCTR in both cases and their amplitude is much larger in the PCTR than in the 
other regions, hence their behaviour is dominated by the PCTR physical conditions. 
The velocity perturbations present a well-defined polarisation, with the amplitude 
of Vx much larger than that of Vz ■ In addition, the perturbations to the density are 
several orders of magnitude larger than those obtained for fast modes (see second 
and third rows of Figure 4). These facts are unequivocal evidence of the slow-like 
character of these modes. In the case Ttr < 7t* the eigenfunctions are very similar 
to those displayed in Figure 9 and have not been plotted here for simplicity. 

As we did for the internal slow modes, one can consider the approximation Vph ~ 
Cstr to the phase speed of the PCTR slow modes. In this way, the dependence of the 

1 /2 

phase speed with the PCTR temperature is given by ~ T^/ . If a similar analysis 
to that of Section 3.3 is performed, one can see that these solutions are very weakly 
affected by the value of the magnetic field strength and the coronal temperature. 
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Figure 9. Perturbations v^, Vx, pi, and Bix (in arbitrary units) versus the dimensionless 
distance to the prominence axis in the case of the PCTR slow fundamental kink mode. The 
top panels correspond to k^Zp = 0.25 (tunnelling branch) and the bottom panels to kxZp = 1 
(propagating branch). The PCTR temperature is Ttr = 500 000 K. The light-shaded region 
corresponds to the prominence slab and the dark-shaded zone to the PCTR. 



5. Discussion and Conclusions 

The investigation performed in tlie present paper is a simple study about the influ- 
ence of the presence of a prominence-corona transition region on tire magnetoacoustic 
osciflatory modes supported by the prominence body. The reader must be aware that 
the wave modes investigated here are not all the possible solutions supported by our 
equilibrium due to the restriction ky = 0. For this reason, the case with ky ^ 
should be investigated in a future work. 

This study is divided into two parts. In the flrst one (Section 3), the JR92 
prominence model for ky = (ER82 case) has been revised and the oscillatory 
modes have been reclassified according to their magnetoacoustic properties. This 
has allowed us to find a wave mode with slow-like properties (the external slow 
mode), which passed unnoticed in ER82 and JR92. In the second part of the paper 
(Section 4), an isothermal PCTR has been considered in the system, and its effect 
on the previously existing wave modes has been assessed. In addition, a new solution 
has been obtained (the PCTR slow mode), which possesses slow-like properties and 
produces large longitudinal motions inside the transition region. 

The main conclusions of this paper are summarised next. 

• Internal slow modes are dominated by prominence conditions whereas the be- 
haviour of the external slow mode is controlled by coronal physical conditions. 
The latter only exists as an evanescent-like solution for small wavenumber. 

• Fast modes are dominated by the prominence but show a slight dependence on 
coronal properties. The influence of the corona on fast modes grows for small 
kx (large wavelength) since perturbations are poorly confined. 

• The presence of an isothermal PCTR does not greatly affect the frequency of 
the oscillatory modes of the equilibrium without transition region if the PCTR 
sound speed is less than the internal Alfven speed. In this situation, a new 
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type of slow-like oscillatory modes (PCTR slow modes) which produce large 
displacements of the plasma inside the PCTR appear in an isolated window 
in the phase speed diagram. These new solutions do not interact with the 
previously existing solutions of the system without PCTR. 
• When an isothermal transition region with a tube speed greater than the in- 
ternal Alfven speed is considered, the PCTR slow modes interact with the 
previously existing fast modes by means of couplings. 

The large longitudinal displacements of the plasma inside the transition region 
produced by the PCTR slow modes could be detected as oscillations in the inten- 
sity or Doppler shift of spectral lines associated with typical PCTR temperatures 
{e.g. Pouget, Bocchialini, and Solomon, 2006). As mentioned in the Section 1, there 
are very few observations of oscillations linked to prominences in transition region 
spectral lines, so one must be cautious about the observational values and their 
theoretical interpretation. Despite these remarks, the addition of the PCTR in the- 
oretical investigations of prominence oscillations is important because it gives a way 
to understand observations obtained with wavelengths corresponding to different 
plasma temperatures. The consideration of an isothermal, homogeneous PCTR is 
a very rough approximation to the real situation but this simple theoretical study 
has allowed us to obtain some interesting results. From our point of view, the most 
important of all these results is that global fast modes of the system can interact 
with PCTR slow modes, which are essentially confined in the narrow PCTR. This 
fact could become of especial relevance if a more realistic, smooth PCTR is assumed. 
Then, a slow continuum would be present instead of individual PCTR slow modes. 
The interaction of the fast mode with this slow continuum could have important 
repercussions on the wave behaviour since a resonance phenomenon could occur. 
The consideration of a smooth PCTR is therefore the next step to be undertaken 
in future investigations, the present isothermal case being a simple first approach to 
the problem. 
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